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We elucidate the basic physical mechanisms responsible for the quantum-classical transition in 
one-dimensional, bounded chaotic systems subject to unconditioned environmental interactions. 
We show that such a transition occurs due to the dual role of noise in regularizing the semiclas- 
sical Wigner function and averaging over fine structures in classical phase space. The results are 
interpreted in the novel context of applying recent advances in the theory of measurement and open 
systems to the semiclassical quantum regime. We use these methods to show how a local semiclassi- 
cal picture is stabilized and can then be approximated by a classical distribution at later times. The 
general results are demonstrated explicitly via high-resolution numerical simulations of the quantum 
master equation for a chaotic Duffing oscillator. 
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I. INTRODUCTION 

Ever since the birth of quantum physics, the bound- 
ary between quantum and classical descriptions of nature 
has been the cause of much controversy and debate. Al- 
though few people now believe in the required existence 
of a "large" classical world in which quantum mechan- 
ics is somehow embedded, even for those that accept the 
primacy of a full quantum description, the identification 
of the actual physical processes that allow a quantum 
dynamical system to be approximated - in some limit - 
by a classical dynamical system often remains less than 
clear-cut. 

Initially, quantum-classical correspondence was 
phrased in the context of understanding how the 
fundamental "subatomic" laws of quantum physics 
could possibly be compatible with a "macroscopic" 
world which, to a very good degree of approximation, 
evolves according to classical Hamiltonian dynamics 
and lacks (classically) bizarre quantum characteristics 
such as interference and entanglement This view 
was famously, if somewhat vaguely, canonized in Bohr's 
Correspondence Principle. The phrase is typically 
invoked to mean one of three related, but not identical, 
subjects: the existence of a formal analogy between 
certain preferred classical dynamical variables and quan- 
tum observables; the limit of large quantum numbers, 
large action or small h, possibly in some combination; 
or the extent to which classical and quantum dynamical 
evolutions agree, in the spirit of Ehrcnfcst's theorem and 
semiclassical dynamics. 

The last two interpretations, which are the principal 
foci of this paper, often overlap with one another but are 
not identical. As an example, the position and momen- 
tum expectation values of a quantum harmonic oscillator 
evolve exactly according to the classical Liouville equa- 
tion, and, given an initial distribution acceptable both 



classically and quantum-mechanically, the two theories 
give identical results. However, when comparing a quan- 
tum energy eigenstate of the oscillator to a classical orbit 
at the same energy, ad hoc reasoning must be utilized to 
eliminate rapid quantum oscillations about the classical 
values, an example of the singular nature of the Ti — > 
limit. That said, interference elimination and dynamical 
agreement are often related insofar as decreasing the size 
of h will usually have the effect of altering the scale of 
quantum interference while simultaneously improving the 
timescale of agreement of classical and quantum expec- 
tation values when they are not already identical (e.g., 
the trivial linear case above). 

It has long been recognized that the problems with 
attaining the classical limit are compounded for nonlin- 
ear systems @. Theoretical analysis and experimental 
observation of chaotic systems over the past forty years 
has made it clear that classical chaos is a real-world 
phenomenon that quantum theory should reproduce to 
within experimental accuracy. Under a unitary quan- 
tum evolution, however, any nonlinear dynamical system 
will eventually fail the conditions of Ehrenfest's theorem. 
Quantum expectation values cannot follow classical pre- 
dictions at long times as quantum mechanics does not 
respect the symplectic dynamical symmetry of classical 
mechanics [3j, |J] . The dynamics of closed bounded quan- 
tum systems are also quasi-periodic; such a system can 
never be chaotic for any non-zero value of Ti. 

A chaotic classical phase space evolution generates 
structures at infinitesimally small scales, whereas, due to 
interference effects, the corresponding quantum evolution 
does not possess a notion of local phase space structures. 
The net effect is a short-time disagreement between semi- 
classical and classical evolutions, followed by a failure of 
the semiclassical approximation itself at longer, but still 
finite timescales ||. This prompted some early in- 
vestigators in the field to wonder if quantum mechanics 
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had to be modified in order to produce chaos Q. As a 
consequence of these obstructions, chaotic systems have 
emerged as a testing ground for whether or not quantum- 
classical correspondence is truly a valid concept, and, if 
so, how it should be properly phrased and addressed. 

A parallel set of experimental developments, particu- 
larly in the last twenty years, have also strongly suggested 
the need for a more refined view of the quantum-classical 
transition (QCT). The border between the macro- world 
of classical mechanics and the micro-world of quantum 
physics has been blurred by technological, observational, 
and theoretical progress. Precision measurements in 
nanomechanics, atomic and molecular optics, and quan- 
tum information processing and communication have 
probed mesoscopic regimes, necessitating a careful anal- 
ysis of the relative merits of using a classical or quantum 
description since the systems studied are neither "very 
large" nor "very small" . In a quite different realm, re- 
cent observations of the cosmic microwave background 
and the large-scale distribution of galaxies have strongly 
supported the notion that primordial quantum fluctua- 
tions seed the formation of large scale structures in the 
Universe @, demonstrating that crude criteria of "mi- 
croscopic" vs. "macroscopic" are no longer sufficient as 
an underlying basis for a serious study of the QCT. Un- 
derstanding the physical mechanisms which define when 
a system behaves classically is now a practical issue. 

A consensus is forming that spanning the gap between 
the above problems and the correspondence principle re- 
quires a robust understanding of open quantum systems 
and quantum measurement [9(. Any experimentally rel- 
evant system is, by definition, a measured system which 
interacts with its environment, if only through a meter. 
A quantum measurement differs from a classical one in 
at least two regards: (i) the intrinsic barrier imposed by 
the uncertainty principle on the precision of phase space 
information a meter can extract and (ii) the more severe 
manner in which the subsystem becomes entangled with 
its environment. Due to this entanglement, quantum 
measurement is generically associated with an irreducible 
disturbance on the observed system (quantum "backac- 
tion"). The desired measurement process must yield a 
limited amount of information in a finite time in order to 
yield dynamical information without strongly influencing 
the dynamics. Hence, simple projective (von Neumann) 
measurements are clearly not appropriate because they 
yield complete information instantaneously via state pro- 
jection. But this fundamental notion of measurement can 
be easily extended to devise schemes that extract infor- 
mation continuously [lol ]. 

The basic idea is to have the system of interest inter- 
act weakly with another (e.g., atom interacting with an 
electromagnetic field) and make projective measurements 
on the auxiliary system (e.g., photon counting). Because 
the interaction is small, the state of the auxiliary system 
gathers little information regarding the system of inter- 
est, and this system, in turn, is only perturbed slightly by 
the measurement backaction. Only a small component of 



the information gathered by the projective measurement 
of the auxiliary system relates to the system of interest, 
and a continuous limit of the measurement process can 
be taken. One then studies the master equation for the 
evolution of the subsystem density matrix conditioned 
on its measurement record. The master equation can be 
further "unraveled" into nonlinear stochastic trajectories 
for a pure state, the so-called quantum trajectories 
An average over the pure states gives back the origi- 
nal density matrix. Unlike in the classical case, where 
the analogous situation refers to a weighted ensemble of 
phase-space points uniquely determined by the probabil- 
ity distribution, a mixed-state density matrix does not 
have a unique decomposition in terms of state vectors. 

It is essential to distinguish between closed evolution, 
where the system state evolves without any coupling to 
the external world, unconditioned open evolution, where 
the system evolves coupled to an external environment 
but where no information regarding the system is ex- 
tracted from the environment, and conditioned open evo- 
lution where such information is extracted. What we 
call the strong form of the QCT describes how a local 
trajectory level picture arises from a conditioned evo- 
lution. However, in many situations, only a statistical 
description is possible even classically, and here we will 
demand only the agreement of quantum and classical dis- 
tributions and the associated dynamical averages. This 
defines the unconditioned weak form of the QCT which is 
the focus of the present paper (for a review see Ref. [13]). 

While the specific nature of the subsystem- 
environment interaction depends on the subsystem 
studied, the actual process of information extraction, 
and unavoidable coupling to other environmental chan- 
nels, there do exist simple, yet physically significant, 
general cases. The systems studied in this paper can 
be interpreted as undergoing a continuous position 
measurement (l3j . where either the results of measure- 
ment are not recorded, or all of the measurements in 
an ensemble are averaged over to erase the information 
regarding specific measurements. Nevertheless, the 
entanglement between the position measuring readout 
and the subsystem still produces a quantum backaction 
in momentum. The form of this open system interaction, 
which falls into the class of Lindblad superoperators, 
rigidly separates the subsystem and its environment [H[ • 

Although a classically chaotic system cannot approxi- 
mate a closed quantum system via the traditional Ti — > 
route, there is good numerical evidence - at least for some 
systems - for the weak form of the QCT. Numerical stud- 
ies of the Duffing oscillator and other systems have shown 
that expectation values of a quantum system subject to 
an unconditioned continuous position measurement will 
come into agreement with the expectation values of an 
(equivalent) open classical system, and that the quan- 
tum phase space will come to capture certain classical 
phase space features [l5j]. In the case of the strong form 
of the QCT, studies have demonstrated the existence of 
nonzero Lyapunov exponents for conditioned systems, as 
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well as inequalities which clearly delineate when the clas- 
sical trajectory interpretation is valid in the conditioned 
case [Hill]. 

An important distinction between the weak and strong 
forms of the QCT must be made. In the conditioned 
case, the master equation actively localizes the wavefunc- 
tion about its expectation value, allowing trajectory level 
agreement between measured classical and quantum sys- 
tems. However, in the unconditioned case, the inequali- 
ties governing the strong classical limit need not be satis- 
fied and localization need not occur. The problem of un- 
derstanding how classical and quantum systems begin to 
look like one another in a generic open system, even with- 
out the advantage of conditioning, has remained open. 

As a final point, we note that while the strong form 
of the QCT must hold for all dynamical systems with a 
classical counterpart, it is not that the weak QCT must 
also do so. The quantum delta-kicked rotor provides a 
particular example of the failure of the weak QCT [l8| . 
The general problem of knowing in advance what gov- 
erns this behavior is not yet resolved, although the work 
in this paper suggests that (effective) compactness of the 
accessible phase space plays an important role. More- 
over, the violation of the conditions necessary to establish 
the strong form of the QCT need not prevent the exis- 
tence of a weak QCT. Since the strong form of the QCT 
requires treating the localized limit, a cumulant expan- 
sion for the distribution function immediately suggests 
itself (l6| . whereas, for the more nonlocal issues relevant 
to the weak form of the QCT, a semiclassical analysis 
turns out to be natural, as will be demonstrated here. 

In this paper we investigate the physical mechanisms 
responsible for the weak quantum-classical transition in 
a one-dimensional, open system with a bounded clas- 
sically chaotic Hamiltonian, expanding on the themes 
of a shorter paper [19(. These arguments are topolog- 
ical in nature and should be generic for compact, one- 
dimensional hyperbolic regions, as well as for unbounded 
systems which stretch and fold in a manner analogous to 
bounded chaotic systems, unlike other studies which fo- 
cus on calculations for a particular system of interest [20| . 
We show how the classical limit is recovered via two 
parallel processes. First, environmental noise modifies 
chaotic classical phase space topology by terminating the 
production of small scale (late-time) structures. (This 
behavior has some parallels with recent numerical studies 
of a chaotic advection-diffusion problem with a periodic 
velocity field, as will be discussed later [Hj]). Second, in 
the quantum picture, environmental noise acts as a reg- 
ulator, attenuating nonlocal contributions to the semi- 
classical wavefunction, and, thereby, stabilizing a local 
semiclassical approximation from the pathologies which 
a classically chaotic system typically generates, so that it 
can now be associated with a noise-modified (smoothed) 
classical phase space geometry. As a consequence of these 
processes, the local semiclassical approximation becomes 
stable at long times, allowing classical and quantum open 
systems to be brought into dynamical agreement at the 



level of distribution functions, rather than the trajectory 
level agreement one obtains from conditioning due to 
measurements. 

The above arguments are very general and apply to a 
wide class of open systems. The key philosophy of our 
approach is that, for a classically chaotic system, corre- 
spondence is inseparable from some notion of measure- 
ment or environmental coupling. We investigate the asso- 
ciated open-system quantum/classical agreement by em- 
ploying the Wigner representation of the quantum den- 
sity matrix and comparing it to the classical phase space 
distribution function, an approach with certain mathe- 
matical and formal advantages We then utilize this 
analysis to elucidate the mechanism by which the agree- 
ment occurs, as well as derive a timescale after which the 
agreement becomes stable. 

Using numerical simulations, we demonstrate the ex- 
istence of the weak QCT for the Duffing oscillator and 
place it in the context of other numerical studies. Many 
of the detailed features of the weak QCT for the Duffing 
oscillator can be explained and predicted by our theo- 
retical framework. We will begin, however, by briefly 
reviewing the semiclassical and classical limits of closed 
nonlinear systems in the Wigner representation, empha- 
sizing why they disagree with their associated classical 
distribution functions at short times and fail as t — > oo. 
For additional background on this topic see Ref . 23] . 



II. THE NONLINEAR CLASSICAL LIMIT IN 
PHASE SPACE 

The Wigner function, /i4/(g,p, t), is a representation of 
the quantum density matrix operator, p, in a c-number 
phase space [2J|. Along with the analogous classical 
phase space distribution function, fc(<l,P,t), we use it 
to compare the dynamics of open quantum and classical 
systems. Using the Wigner function as a tool for study- 
ing the quantum-classical transition is conceptually and 
practically advantageous. It allows one to compare classi- 
cal and quantum dynamics in phase space (though there 
are pitfalls one must be aware of), rather than trying 
to compare, say, wavefunctions in L 2 to classical trajec- 
tories. More importantly for our purposes, the theory 
of semiclassical approximations can be directly tied to 
the evolution of classical curves in phase space, making 
it easier to visualize the extent to which quantum and 
classical dynamical evolutions agree [22]. For a classi- 
cally chaotic system, distribution functions can also give 
a clearer sense of global phase space topology, allowing 
one to examine the extent to which dynamical agreement 
over an entire compact hyperbolic region of interest is 
achieved. 

The Wigner representation of an operator, A, is defined 
as: 

/oo XX 
dXe-^ x ' n {q+-\A\q--). (1) 
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The Wigner function is the Wigner representation of the 
general mixed state density operator p = Ci\ipi) (ipi\ 
yielding 

fw(q,P,t) = 7^ / dX exp ( - j x 



2ttTl ./_™ \ Ti 

It follows that, for any operator. 



J2ciiPi(q+^, m(q-?-,t). (2) 



(A) = Tr(ip) 

/ dp A w (q,p,t)f w (q,p,t). (3) 
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Any classical quantity, q n p m , can be associated with a 
quantum operator, via the Weyl ordering 



1 oo 



(4) 



r=0 



Thus one can compute averages of any classical quantity 
in the Wigner picture. 

Unlike a classical phase space distribution function, the 
Wigner function is only a quasiprobability distribution, 
as it can take on negative values. This condition also im- 
plies that the Wigner function cannot generally be used 
as a conditional probability distribution and is bounded 
by ±(7r7i) , which prevents it from being a delta func- 
tion in phase space at finite h, and, therefore, prevents 
it from representing a classical trajectory (25[. The clas- 
sical phase space distribution is a true positive definite 
probability distribution capable of determining averages 
over arbitrarily small phase space regions, whereas the 
degree to which a Wigner function can capture a local 
average depends on whether the region being integrated 
over is receiving strong quantum interference effects from 
locations outside of the integrated region. 

The equation of motion for the Wigner function is 
given by the Wigner representation of the equation of 
motion for the density operator: 



at 



— L c fw + LqJv, 



where the classical Liouville operator 
L c 

and the quantum correction 



dV 
pd x + dp 
ox 



L 



Q 
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J (2 2 ™(2n + l)!) 



d 2 x n+1 Vd 2 p n+1 . 



(5) 



(6) 



(7) 



The form of this evolution equation suggests an intuitive, 
but misleading, interpretation of the how the classical 
limit is achieved [26| . In the equation of motion, h only 
appears in the L q term. So it is tempting to suggest that, 



as h — > 0, the "quantum contributions" to the evolution 
of the Wigner function likewise decrease. However, all of 
the momentum derivative terms in both Lq and Lc are 
proportional to h 2n d 2n+1 f w . Since, by definition, 



f w ~ exp 



(8) 



after one takes the appropriate momentum derivatives, 
it is clear that, like the wavefunction, fw is 0{fi^ 1 ) to 
leading order in Ti. One can never expect quantum correc- 
tions to smoothly disappear as h is decreased due to this 
essential singularity, which produces increasingly rapid 
oscillations as ti — > 0, and will keep the Wigner function 
from tending to a positive distribution. To eliminate the 
rapid oscillations, one often introduces an ad hoc filter, as 
in the case for the Husimi-type Gaussian filters [13, . 
This can forcibly produce positive-definite distributions 
but lacks an underlying dynamical justification. 

The formal study of the classical limit in phase space 
begins by constructing a semiclassical Wigner function 
from an underlying semiclassical wavefunction. The 
semiclassical wavefunction is the singular and 
constant part of a general wavefunction in the Ti — * 
limit [1^]. In this sense, any small h view of the classi- 
cal limit must focus on the semiclassical regime since the 
semiclassical wavefunction is the irreducible part of the 
wavefunction in this limit. The standard presentation 
tends to view this process as simply representing the two 
lowest order terms in a perturbation series for the phase 
of the wavefunction. However, the higher order terms 
in this series, in addition to being notoriously difficult 
to calculate, arc rarely useful. The remaining terms can 
be thought of as a vanishing, 0(h) error, and not as a 
series of hig her order terms waiting to be explicitly cal- 
culated [3Cj |. 

Most importantly, a semiclassical wavefunction is di- 
rectly associated with the evolution of classical phase 
space curves. The formal procedure constructs an initial 
wavefunction from an TV-dimensional Lagrangian mani- 
fold embedded in a 2 TV-dimensional phase space [3(| ■ In 
this paper, phase space is two-dimensional and so the as- 
sociated Lagrangian manifold studied is a curve, which 
will be one of a family of phase space curves parametrized 
by the continuous parameter 7, as elucidated in Ref. [Hj]. 
An initial semiclassical wavefunction associated with the 
curve 7 will have the form: 



^(<2,0;7) = A) (<7; 7) exp 



-So(q;i) 



(9) 



where Ao(q; 7) and So(q', 7) are real- valued functions. For 
simplicity of presentation, we will assume that this ini- 
tial curve has a single momentum value associated with 
each position. Relaxing this assumption would result in 
a slightly more awkward presentation, but would not al- 
ter its substance. The above form naturally induces a 
Lagrangian curve in phase space if the associated mo- 
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mentum has a well-defined classical limit. Namely, 

limaJ-^(<7,0,7) = A (g, 7 ) ^ g ' 7) = p(q,j). (10) 
n—>0 Oq Oq 

The evolving classical curve will typically develop turn- 
ing points, which can result in multiple momentum values 
for a given position. This will certainly be the case for 
the highly nonlinear systems addressed here. As a con- 
sequence of this folding, one assigns a new action to each 
branch of the curve as it evolves, as described in Ref. [22| . 
The action at time t for the j-th path is then given as 

s j(q,t\l) = S (q 0j ;j)+ ( dq'pj(q',t;j) 

Jqoj 

ft 

dt'H(q 0j , p 0j (q j , t'; 7) , t') , (11) 

10 

which yields the semiciassicai wavefunction: 

N , . . 

ip(q,t;~f) = ^2A j (q,t)eKpi-S j (q,t;j)-—, 



i=i 

+o(h), 



(12) 



where Vj is the j-th Morse index, defined as the number 
of times the determinant is equal to zero along the path 
connecting {q oj ,poj) to (q,Pj). 

By substituting the semiciassicai wavefunction into the 
definition of the Wigner function, one can construct a 
geometric interpretation of the accuracy of a semiciassicai 
analysis. For the purpose of clarity, we will assume we are 
dealing with a pure state density matrix, the extension 
to mixed states being straightforward. The semiciassicai 
Wigner function becomes: 

1 f°° 

fw{q,p,t;i) = ^ J dx - (?, t ; 7) x 



exp 



- (q + X/2, t; 7) - S s (q - X/2, t; 7) - pX) 



ITT 



(13) 



where Aij = A t (q + X/2, t; j)A j (q - X/2, t; 7). To get a 
sense of the primary contributions to this integral as h is 
brought to zero and the integrand rapidly oscillates, we 
examine the stationary phase condition: 



— (S z (q + X/2,t) 



S j {q-X/2,t)-pX)=0, (14) 



where, for clarity, the 7 parameter is suppressed for the 
remainder of the paper. In the stationary phase approx- 
imation, the Wigner function is separated into a singu- 
lar stationary part and an additional O^ 1 ^ 2 ) oscillatory 
part [3l|. Therefore, the stationary phases are the most 
relevant contributions in the h —> limit, as rapid os- 
cillations become less significant. After substituting the 



expression for the evolved action, the stationary phase 
condition becomes: 

l -{ Pl {q + X/2,t) + Pj {q - X/2,t)) = p(q,t). (15) 

If i = j, this is the famous Berry midpoint rule: \{pi{q + 
X/2)+ Pj {q-X/2)) = p{q,t) '22\. That is, the stationary 
phase contributions at a point {q,p) come from the aver- 
age of the momenta on a given solution curve evaluated 
at the end of an interval of width X about q. 

If the point (q,p) is particularly close to a curve, Pi(q), 
as will often be the case when the underlying curve 
evolves in a chaotic region of phase space, then the sta- 
tionary phase points will coalesce, invalidating the sta- 
tionary phase method. Likewise, the WKB wavefunction 
itself is not valid near turning points, as the Jacobian van- 
ishes. However, these cases are remedied by the uniform 
approximation, which yields a symmetric Airy function, 
rather than sinusoidal, behavior 32]. Therefore, if (q,p) 
is too close to a given branch, the expression for the semi- 
classical contribution for that branch should be replaced 
by the uniformized form. The only remaining problem 
which can invalidate the expression is the appearance 
of catastrophes when (q,p) is a focal point of a curve, 
which can be dealt with analytically, as was also studied 
by Berry, but is outside the realm of this paper. Close to 
the classical curve Pi(q), the uniformized WKB approx- 
imation for the Wigner function has an Airy "head" of 
width ~ h ' and peak height ~ h~ ' . In this limiting 
case, the uniform approximation can be further simpli- 
fied and written in the form of a "transitional approxi- 
mation" which is valid only very near Pi(q). Remarkably, 
the h — ■+ limit of the transitional approximation is in- 
deed a classical delta function, which allows the Wigner 
function picture to give a physically clearer presentation 
of the classical limit. 

The semiciassicai quantum Wigner function goes 
through three phases in its evolution, if the underlying 
classical dynamics is chaotic, as laid out by Heller and 
Tomsovic ja]. During the (very short) first phase, if the 
initial condition used is that of a classical distribution, 
there will be little disagreement between the quantum 
and classical evolutions. This is followed by a second 
phase, where the semiciassicai approximation reproduces 
the wave function dynamics, but is distinctly nonclassi- 
cal. At a longer timescale, proportional to inverse powers 
of h, the semiciassicai approximation fails, as the dis- 
tance between classical manifolds becomes so close that 
the cumulative interference cannot be locally ascribed to 
any given curve. So, in the first, classical regime, there 
is little interference. In the second, semiciassicai regime 
there is some, possibly strong, quantum interference, but 
it is in the form of local fringing about classical curves. 
In the final, fully quantum phase, there is strong global 
interference, and local classical manifold evolution is of 
little relevance to the quantum propagation. 

There are two sources of quantum interference in phase 
space: local Airy "shadows" of the short wave classical 
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curve and nonlocal contributions from multiple curves, 
the latter being more problematic if we wish a weak QCT 
to hold. In order to maintain a stable classical limit for 
a classically chaotic system, it must be possible to keep 
a system in a more or less classical regime (analogous to, 
but not the same as the first regime discussed above), 
allowing only a small admixture of local interference ef- 
fects. We show below how such a stable classical limit 
arises in open systems, via the same physical process that 
simultaneously leads to a smoothing of the classical phase 
space geometry. 



III. OPEN SYSTEMS AND MEASUREMENT 

To model the interaction between a subsystem and its 
measuring device we choose the form of an unconditioned 
continuous position measurement. This provides the 
minimum level of interaction necessary to bring quantum 
and classically chaotic dynamical systems into (approxi- 
mate) agreement with one another at the level of distri- 
bution functions. The model of a conditioned continu- 
ous position measurement (i.e., evolution of the system 
density matrix taking the results of measurement into 
account) is given by the following master equation [33| : 



dp = --[H,p] + k[X[X,p]]dt 

-^([X,p] + -2p(X))dW, (16) 
where the observed measurement record is given by 



evolution: 



dy = (X)dt+ T dW. 



(17) 



In the above equation (X) = Tr(pX), dW is the Wiener 
measure [(dVF) 2 = dt], k represents the strength of the 
interaction between the subsystem and the measuring ap- 
paratus and k measures the rate at which information 
about the system is being extracted. The fractional mea- 
sure of extracted information is given by the efficiency of 
the measurement r\ = k/8k. The first term in Eqn. (|16|) 
is just the unitary evolution for the closed system, the 
second is a diffusive term arising from quantum backac- 
tion, and the third represents the conditioning due to the 
measurement. 

The conditioned evolution can localize the state about 
the measured position value; the extent of this localiza- 
tion (proportional to k) must however be tempered by 
the associated increase of backaction noise (concomitant 
increase in k). Nevertheless, inequalities can be derived 
that show under what conditions both of these conflicting 
effects can be reconciled and agreement between classical 
and quantum dynamics achieved at the level of trajecto- 
ries [TH, [TtJ - the strong form of the QCT. 

If one averages over all obtained measurement records, 
one obtains the master equation for an unconditioned 



dp=--[H,p] 



k[X[X, p]]dt. 



(18) 



This evolution can also be achieved by setting the effi- 
ciency of the measurement, and, therefore, k = 0. Once 
one does so, the localization inequalities which character- 
ize the strong form of the QCT fail, showing the inability 
of the weak QCT to capture trajectory level chaos and 
the need for the distribution function approach employed 
here. The evolution equation is the same as that for 
the Caldeira-Leggett model in the weak coupling, high 
temperature approximation [34]. The key point here is 
that while the conditioning term is absent, the backac- 
tion term remains. This is very different from the classi- 
cal case, where averaging over measurements simply gives 
back the closed-system Liouville equation, thus highlight- 
ing the contrast between the active nature of quantum 
measurements versus the passive nature of classical mea- 
surements. 

The master equation (fT5| is the starting point in our 
analysis of the weak QCT utilizing the Wigner function. 
In the Wigner representation, this equation becomes 



Of 



ir 



dt 



Lcfw + Lnfi. 



QJW 



D 



dp 2 



\v 



(19) 



where the diffusion coefficient D = Ti k. If we set Lq — 
0, we obtain a dual classical evolution equation, for the 
classical distribution function fc(q,P,t): 



dfc 
dt 



~et — Lcfc 



D 



d 2 f c 
dp 2 



(20) 



which, given its form, we will call the dual Fokker-Planck 
equation. Note that this Fokker-Planck equation does 
not represent the dynamics of an associated classical ob- 
served system. Here it has two key roles: it represents the 
classical template for a semiclassical open-system anal- 
ysis and also the proper (approximate) classical limit- 
ing form if the weak QCT were to hold. This particular 
Fokker-Planck equation is better viewed as simply a clas- 
sical dual of the quantum master equation (|19|) , without 
an independent physical existence. 

A final note on timescale separations is necessary to 
clarify the physical situations under which Eqns. (|18p and 
(UHJ) are considered to hold. We are not interested in im- 
posing initial conditions on the quantum dynamics that 
have classical analogs (e.g., Gaussian wavepackets), and 
then looking for the emergence of short-time quantum ef- 
fects. In fact, we acknowledge the existence of quantum 
initial conditions explicitly (as in the numerical simula- 
tions of Section VI), and investigate quantum-classical 
convergence in the sense of the convergence of distri- 
bution functions as obtained from the quantum master 
equation and its classical dual. 

At the the same time, we are particularly interested 
in the dynamics set by the closed-system Hamiltonian, 
with minimal influence from the external environment 
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or continuous measuring process, i.e., the weak coupling 
limit. In this limit, we can ignore the dissipative ef- 
fects of external couplings (clamping due to environment 
modes and/or measurement backaction), but consider 
only diffusive effects, which remain finite in the weak cou- 
pling limit (as in the weak-coupling, high-temperature 
Caldeira-Leggett model). There are two timescales as- 
sociated with these statements. The first, t re i axi is the 
time taken for the system to relax to a thermal state, 
(or noncquilibrium steady state depending on the cir- 
cumstances) and is typically controlled by the matching 
of energy exchange as set by the dissipation and diffu- 
sive channels. The second, tdiff, is the diffusive heating 
timescale which, in the absence of dissipation, leads to 
continuous heating of the system. Since at late times, 
when dissipative effects would be expected to occur, this 
heating is unphysical, it is clear that our analysis assumes 
t <C tdiff- Therefore, we are interested in the dynam- 
ics of open quantum systems on intermediate timescales, 
longer than the system dynamical timescales, yet far from 
the (asymptotic) timescales relevant for close to steady- 
state behavior. All remarks below on "long-time" behav- 
ior apply to this intermediate timescale and not to some 
eventual steady state. 



IV. MODIFICATION OF PHASE SPACE 
GEOMETRY FOR A CHAOTIC SUBSYSTEM 

The first step in our analysis is the study of the dual 
Fokker-Planck equation. As earlier mentioned, following 
a semiclassical line of reasoning, the motivation for this 
is that the measurement / environmental interaction mod- 
ifies the geometry of a chaotic classical phase space in a 
manner which can allow dynamical agreement between 
classical and quantum systems. The key point is that, 
due to the diffusion term, one necessarily sees a termina- 
tion in the level at which one can discern the long-time 
development of fine structure. The (exponential) long- 
time development of structure is a hallmark of classically 
chaotic systems in a compact space, and, as discussed in 
the previous section, leads to disagreement between clas- 
sical and semiclassical results, followed by a complete 
failure of semiclassical analysis. But, as these structures 
are averaged over, the resulting smoother phase space 
geometry can be consistent with the existence of a local 
semiclassical description. 

We will show below that the diffusion term in the 
Fokker-Planck equation terminates the development of 
small scale structures at a finite time, denoted by t* . 
At this time, there will be an associated area, l c i(t*) 2 , 
below which no smaller phase space structures can be 
discerned. To understand the termination of structure, 
we consider the Langevin equations underlying the dual 
Fokker-Planck equation. These are given by 

dq=pdt/m (21) 



and 

dp = f{q)dt + V2DdW, (22) 

where f(q) = —dV(q)/dq, dW is the Wiener measure 
[(dW) 2 — dt], and D is the noise strength. Since D is 
constant, one can consistently write dW = £(t)dt, where 
is a rapidly fluctuating force satisfying = 

and (£(£)£(£')) = 6(t — t') over noise averages. 

A hyperbolic region of the phase space of a bounded 
chaotic Hamiltonian system is foliated by its unstable 
manifold, which emerges from the stretching and folding 
behavior induced when nonperiodic solution curves are 
confined to a bounded region. A trajectory in the neigh- 
borhood of a hyperbolic fixed point will create large scale 
structures, due to its exponential growth away from the 
hyperbolic point. As it evolves, since it can only explore 
the energetically allowed region, it will fold onto itself and 
create smaller scale structures. For a bounded chaotic 
region, the curve will eventually fill the allowed space. 
The important consequence for this analysis is that this 
filling is done preferentially. Large scale structures are 
initially generated by rapid stretching and are associated 
with short timescales. The smaller scale fine structures 
are then filled in afterwards as the system continues to 
fold on itself and are, therefore, a late-time feature. 

In order to investigate how environmental noise mod- 
ifies this picture, we perform a perturbative expansion 
of the solution curve in the small noise limit in the 
neighborhood of a hyperbolic fixed point (q eq ,0), where 
f(leq) = 0, where \/2D is treated as the small noise 
parameter [35L [36[ . As emphasized in the previous sec- 
tion, this assumption is physically justified by the argu- 
ment that the affected noise scale in phase space should 
be smaller than that of the system dynamics. As dis- 
cussed in the introduction, the system is taken to be 
weakly interacting with its environment to ensure that 
the system dynamics arc affected only pcrturbativcly. 
To leading order in V2D, we can therefore separate the 
dominant systematic components from the noisy compo- 
nents via q(t) sa q c (t) + qff(t) and p(t) w pc(t) +PN&), 
leading to the usual Hamilton's equations for qc and 
Pc, and to the coupled equations dq^ = p^dt/m and 
dpN = m\ 2 qjydt + dW, where m\ 2 = df(q eq )/dq defines 
the local Lyapunov exponent, A. These have the solution, 

q(t) = q eq +C+e xt + C-er xt 

with an analogous expression for p(t) . 

To understand the effect of noise on the foliation of 
the unstable manifold, one needs to transform from the 
position and momentum basis into the stable and unsta- 
ble directions. The dimensional scalings q' = \J Xmq and 
p' = pj \J Xm are introduced so that the rescaled posi- 
tion and momentum have the same dimensions and also 
so that the stable and unstable directions are orthogo- 
nal. An arbitrary time rescaling, which would give the 
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correct units, would not guaran tee o rthogonality. If we 
project the solutions for q' = \J Xmq and p' — pj\JXm 
along the stable (-) and unstable (+) directions, we find 
the following expression for the components of the noisy 
trajectories evolution in these two directions: 



u±(t) = 4§(<?'±f>') 

= V2W7±e ±At ± 




du£(u)e ±A (*-^4) 
(25) 



One can now analyze the effects of these noisy trajec- 
tories on the evolution of the distribution function which 
they unravel. The average over all noisy realizations of 
the displacement in the stable and unstable directions is 
given by (u±) = \ / 2XmC±e ±xt , as expected from a per- 
turbation in the neighborhood of a hyperbolic fixed point. 
More information is found in the second order cumulants. 
Whereas, the stable and unstable directions have vari- 
ances of ±(D / (2mX 2 ))(e ±2Xt — 1), the off-diagonal cu- 
mulant is (it + ii_) — = —Dt/(mX), displaying 

the linear spreading associated with a Wiener process. In 
forward time, where the evolution of a trajectory is de- 
termined by the unfolding of the unstable manifold, this 
spreading indicates that, as the trajectory evolves, it will 
simultaneously smooth over a transverse width in phase 
space of size 



l d (t) « y/Dt/(m\). 



(26) 



One is left with a picture of a curve following a classi- 
cal path in the unstable direction while carrying small 
amounts of transverse noise. In a bounded, compact 
phase space region, this implies a termination in one's 
ability to measure the position and momentum of the 
trajectory on a scale smaller than the aforementioned 
width. In other words, the fine structures associated 
with a chaotic region will be smoothed over in the av- 
eraging process, causing the development of large scale 
structures which occur prior to this termination time to 
become pronounced. 

Given a set of parameters associated with this compact 
phase space region, one can estimate the value and scal- 
ing associated with the termination time, t*. Consider 
an initially small compact region of phase space area Uq, 
then its current phase space "length" is approximately 
uoe xt , where A is the time-averaged positive Lyapunov 
exponent. If the trajectory is bounded within a phase 
space area A, the typical distance between neighboring 
folds of the trajectory is estimated by 



6(t) » A/(u e xt ). 



(27) 



This formula only applies once the curve has begun to 
fold on itself; given a particular choice of uq, one must be 
careful that enough time for folding to occur has passed 
before using the above equation. One can, in this spirit, 



estimate a "folding time" and compare it with the even- 
tual computed value of i* to again insure that this anal- 
ysis is self-consistent. In any case, one cannot make Uq 
arbitrarily small when exploring the QCT because the 
uncertainty principle sets a lower bound on phase space 
area. Note that the length of a long timescale - long com- 
pared to the dynamical timescale -is also implied by the 
appearance of the time-averaged Lyapunov exponent, A. 

Phase space structures can only be known to within 
the width specified by the noisy dynamics, hence there 
will come a time at which the rapidly falling scale 6(t) set 
by the folding will be smaller than the slowly increasing 
filter scale, l c i(t), at which lengths are averaged over as 
given by Eqn. (|26|) . The time at which any new structures 
will be smoothed over is given by equating Eqns. (|26[) and 
(j27)l . to yield 



Dt* 



A 



(mA) Uq 



exp(— Xt*), 



(28) 



a simple transcendental equation for t*. Typically, t* 
is expected to be significantly larger than 1/A, as there 
will usually be many foldings before the filtering becomes 
effective. In this case a simple iterative procedure can be 
used to find the approximate solution, 



t* 



2A 



1 



ln(x ) 



(29) 



where x = ln(2mA 2 A 2 / '(DuD) . 

After the time t* no new structures will be discerned, 
since they will be smaller than the averaging scale set by 
the noisy dynamics. This implies the existence of a phase 
space area l c i(t*) 2 below which phase space structures are 
smoothed over. As a result, the dual Fokker-Planck equa- 
tion for a chaotic system is such that we can only discern 
large scale structures (small and large being relative to 
the cutoff l c i(t) produced prior to t*). When construct- 
ing a classical limit for the open quantum evolution (|19|1 . 
we now only have to capture the larger, short-time dy- 
namical features and not the full chaotic evolution of the 
classical Liouville equation with its - from a quantum 
perspective - small-scale pathologies. 

Before proceeding further, we mention an analogous 
situation in studies of chaotic advection-diffusion in fluid 
dynamics. The evolution equation for the concentration 
density, c(x, t) , of a set of particles diffusing in a fluid 
without sinks or sources is given by 



dc 

— + Vc • v = kV 2 c, 
at 



(30) 



where v(x, t) is the velocity field of the tracer particles. 
This matches the classical Fokker-Planck equation stud- 
ied here, if one sets 



P_ 

rn ' 



dV 
dq 



(31) 



and the gradient is taken with respect to q and p. Dif- 
fusion, in our case, is only with respect to p. The phase 
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space distribution function is then regarded as the con- 
centration of particles in phase space in a given region, 
which is certainly an appropriate interpretation. A nu- 
merical analysis performed in Ref. [21] showed results 
similar to our predictions where, for a certain vale of k, 
equivalent to D in our case, the evolution converged to 
a stationary pattern at a finite time, with only residual 
diffusion afterwards. The final pattern was termed an 
inertial manifold and related to the unstable manifold. 
Following this, the existence of such a manifold beyond 
a critical k value was demonstrated analytically [371 ]. 

This suggests that the qualitative classical analysis 
provided here might be made more rigorous. The anal- 
ysis of Ref. [13] relied significantly on applying peri- 
odic boundary conditions to the concentration evolu- 
tion and exploiting the resulting gaps in the spectrum 
of the Laplacian. The open boundary conditions rele- 
vant to this paper, however, appear to preclude such an 
approach. Nevertheless, qualitative similarities do exist 
and the two fields may well inform each other in future. 
(Of course, this is a purely classical analysis and does not 
bear directly on the quantum evolution, except implicitly 
since semiclassical evolution tracks the classical manifold 
structure.) 



V. SEMICLASSICAL ANALYSIS FOR AN OPEN 
CHAOTIC SYSTEM 

We now turn to the semiclassical analysis of the open 
system master equation (TT9")) in order to estimate the 
conditions under which a weak QCT might exist. We 
begin by rewriting the semiclassical Wigner function in 
the weak noise limit utilized in the previous section. In 
this limit, the classical action is modified to S(q,t) « 
S(q c ,t) - V2D f*dt£(t)qc(t), as in Ref. 38]. The first 
term will evolve classically, as discussed in Section III, as 
will the position coordinate which appears in the second 
term. If we insert the above semiclassical action into the 
expression for the Wigner function we get the following 
result: 



classical Wigner function is obtained: 



fw{q,p,t) 



dX 



-ipX/H N 



5> 



2nh 13 



exp (-~XV2D^ dt' £(*'; 



exp ( ~(Sa(qc+,t,P) - S Cj (qc-,t,P)) 



ITT 



(32) 



noting that, since the amplitude is a second derivative 
and the noisy perturbation is linear, noise only effects 
the action to lowest order. 

If we next average over all noisy realizations, the fol- 
lowing suggestive expression for the noise averaged semi- 



dX 



^—ipX/h 

2irh 



-p(-f^ 2 )f> 



exp ( -{S Ci (qc+,t,P) - S Cj (qc-,t,P)) 



(33) 



The only alteration to the expression for the semiclassi- 
cal wavcfuntion to lowest order in the noise strength is 
the appearance of a new Gaussian term. The presence 
of noise acts as a dynamical low-pass Gaussian filter of 
semiclassical phases, attenuating large X contributions. 
For any solutions to the above equation, phases will be 
suppressed which have wavelengths greater than 



x w h/Voi. 



(34) 



These are the long, nonlocal "De Broglie" wavelength 
contributions to the semiclassical integral, the very sort 
of contributions previously identified as being particu- 
larly problematic in terms of obtaining a weak QCT. The 
filter prevents the integral from becoming overwhelmed 
by long range contributions as stretching and folding oc- 
curs which can lead to disagreement with classical results, 
as well as the eventual failure of the approximation. 

We now combine the above with the classical result 
from the last section. It is seen that the diffusion 
causes two effects: suppression of nonlocal phases in the 
semiclassical integral beyond a certain scale given by 
Eqn. (j3"4"|) and a smoothing of the dual classical phase 
space over fine structures smaller then a scale given by 
Eqn. (|28p . Each of these effects overcomes the two semi- 
classically identified difficulties associated with a weak 
QCT for chaotic systems: the Wigner function is no 
longer dominated by nonlocal contributions and also does 
not need to track, nor does it receive interference from, 
very fine scale structures. From these two scales we 
should, therefore, be able to set a (semiclassical) crite- 
ria for the existence of a weak QCT for a bounded one- 
dimensional chaotic system. Physically, the local semi- 
classical approximation is valid when the primary contri- 
butions to the semiclassical integral at a given point (q,p) 
come from the local branch of the trajectory on which 
the point is located. This will occur only when the scale 
at which local classical smoothing occurs matches or ex- 
ceeds the filtering scale for semiclassical phases. When 
this occurs the nearest possible branch which is capable 
of delivering nonlocal interference effects will have those 
effects filtered within the semiclassical integral. As a re- 
sult one can recover the usual short wave semiclassical 
picture of a trajectory "decorated" only by local interfer- 
ence fringes. 

More specifically, if we rescale the filtering condi- 
tion (|34p in phase space units 



In 



(35) 
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then the semiclassical criterion for the weak QCT is given 
by l q (t) < l c i{t). The quantum scale l q decreases with 
time, as the noise filtering, beginning with "fast" phase 
space oscillations (interference due to far-separated fea- 
tures in phase space), reaches down to ever "slower" 
interference scales (due to small-scale phase space fea- 
tures). If small-scale classical structures continued to be 
generated in phase space at a rate outstripping the de- 
crease of l q with time, the QCT would not occur. Due to 
the presence of noise, however, classical small-scale struc- 
ture does not grow exponentially, but is eventually cut 
off by l c i, which grows with time. Therefore, l q and l c i 
must cross each other, and this point defines the weak 
QCT timcscalc, t qc . For t > t qc , the classical structures 
are large enough that the noise filtering is effective in 
smoothing over the associated interference terms. 

The weak QCT timescale follows from equating 
Eqns. (f!?rj]) and (f3"o]) for l c i and l q , respectively, 

t qc mhX/D. (36) 

Using Eqn. (|26|) once again, one finds that this condition 
is nothing but Q(£ gc ) — l q (tqc) ~ h, which would have 
been suggested by basic intuition. Following from the dis- 
cussion above, t qc can also be interpreted as the timescale 
beyond which a semiclassical approximation becomes sta- 
ble for an open quantum system. After this time, clas- 
sical dynamics should approximate quantum dynamics 
sufficiently. 

Note that the two timescales discussed so far, t* and 
t q ci scale very differently with the diffusion coefficient, 
D. Whereas, t* ~ ln(l/D), t qc ~ 1/D, implying that the 
timescales are far-separated in the small D (weak noise) 
limit, where typically, t qc t*. It is possible, however, 
to have t qc < t* even at modest values of D. The phys- 
ical interpretation of these two possible situations is as 
follows. As discussed previously, the timescale t* sets the 
"freeze-out" of classical phase space structures, but it is 
possible to have a weak QCT occur on either side of the 
freeze-out. When t qc > t*, even though the large-scale 
classical phase space template is relatively fixed, small- 
scale discrepancies will exist between the quantum and 
classical distributions at least until t ~ t qc . Though, in 
this case, the classical filtering will have terminated the 
development of classical structures at t* , some time must 
still elapse before interference between branches has been 
sufficiently filtered. On the other hand, when t qc < t*, 
the weak QCT can occur while the large-scale classical 
phase space structures are still evolving since the classical 
freeze-out has not yet taken place. 



VI. NUMERICAL SIMULATIONS 

The analysis in the preceding sections has helped to 
establish a set of criteria which, once met, allow the ex- 
istence of a weak QCT for classically chaotic systems. 
Given their somewhat heuristic nature, it is important 



to examine these predictions numerically. In the quan- 
tum evolution, once the inequalities are satisfied, noise 
will filter nonlocal quantum interference between the sur- 
viving large scale phase space structures, so large scale 
coherences should not be present. If not, one will es- 
sentially see a global phase space diffraction pattern, 
with large-scale coherences persisting between all parts 
of the bounded phase space region, an example of which 
is shown in Figure 1. 
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FIG. 1: Phase space rendering of the Wigner function for the 
Duffing system at time t = 314 periods of driving. The non- 
local interference is significant and cannot be associated with 
specific classical structures. This plot is taken at a relatively 
small D value (10 -4 ) for resolution purposes. The value of Ti 
is set equal to 1 in order to clearly demonstrate this effect. 



The most direct numerical test is a close examination 
of the time evolution of both the classical and quantum 
distribution functions for the quantum and dual classical 
evolutions. In this manner one can examine whether, 
over the expected timescales as predicted in the previous 
section, the expected phase space features are present for 
the template classical distributions and quantum Wigner 
functions. 

A direct examination is necessary as other, seemingly 
logical measures, can sometimes be misleading. For in- 
stance, looking at expectation values is not always help- 
ful. Typically, one would sample a set of lower order mo- 
ments, and follow their expectation values for a desired 
amount of time. However, any classically chaotic system 
will have, over time, an infinite set of non- vanishing mo- 
ments, as will any non-Gaussian quantum system [3^ |. 
That said, one might attempt to argue that the effect 
of these higher-order moments may well be negligible. 
Even if this were the case, however, numerical simula- 
tions of chaotic systems have failed to find well-defined 
break times at which even lower-order moments perma- 
nently separate, and expectation values can agree well for 
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surprisingly long timescales and even without the pres- 
ence of environmental noise @. Other measures, such 
as suppression of the integrated negativity of the Wigner 
function, also are not necessarily signatures of quantum- 
classical correspondence, as shown in Ref. |4|. Once the 
amount of negativity is eliminated via environmental de- 
coherence, it was shown by these authors that the quan- 
tum evolution of a system may still disagree with its clas- 
sical counterpart. 



A. Numerical Methods 

Numerical solutions of the quantum master equation 
(TTU|) for the Wigner function and of the corresponding 
dual classical Fokker-Planck equation were carried out 
using a split operator spectral method implemented on 
parallel supercomputers [4(|. The spectral method is 
particularly well-suited to high spatial resolution simula- 
tions where spatial structure is cut off above some given 
wavenumber - this is the case here for both the quantum 
and dual classical evolutions. 

The time-stepping strategy is the same as that in anal- 
ogous classical symplectic integrators. Suppose the time 
evolution of a function, f(t), satisfies the operator equa- 
tion: 



^l = (L A + L B )f, 



(37) 



where the separate evolutions given by La and Lb can be 
implemented exactly. The exact solution to this equation 
is given by: 



/(*) 



(L A +L B )t 



/(o). 



(38) 



momentum diffusion operator. As each piece involves ei- 
ther derivatives of position or momentum, but not both, 
the individual operators can be easily evaluated using a 
fast Fourier transform. The split-operator method pre- 
serves the unitarity of evolutions when D = and given 
a sufficient number of grid points in the spatial and mo- 
mentum directions - satisfying associated Nyquist condi- 
tions - the operators can be evaluated at each timestep 
with essentially no spatial discretization error. 

The typical mesh used over phase space consisted of 
4096 by 4096 grid points. This size was determined by 
our need to resolve the bounded classical phase space por- 
trait for the amount of time necessary to show long range 
agreement between the classical and quantum evolutions. 
If D = 0, the classical phase space will be chaotic, and 
the system can only be explored for short times, before 
which structures begin to proliferate on scales smaller 
than the area defined by the grid spacing. The addition 
of an environmental interaction, as demonstrated in the 
theoretical section, prevents structures from forming on 
infinitely small scales. This makes it possible for the clas- 
sical evolution to converge as resolution improves. The 
aforementioned grid-size is the one for which convergence 
was achieved for the systems we studied, and was derived 
empirically. Convergence for the quantum evolution is 
determined by the smallest scales - Sx = Ti/P in space 
and dp — h/L in momentum - present in the Wigner 
function (L and P are the scales of the system boundaries 
in length and momentum, respectively). Thus, a typical 
mesh spacing can be fixed without regard to the strength 
of the environmental interaction. For our investigations, 
this required less resolution than in the classical phase 
space and, therefore, the dual classical evolution dictated 
the grid-size for the numerical simulations. 



Since La and do not commute in general, the fact 
that the individual evolutions are known exactly is not of 
direct use. An integration scheme for a small timestep At 
can be constructed simply, however, using the Campbell- 
Baker-Hausdorff theorem: 



/(At) 



(^L A ) p (AtL B ) e (^L A ) 



e \ 2 



>/(0) + O(Ai 3 ). (39) 



With the assumption that the exponentiated operators 
can be applied exactly, this method is accurate to second 
order in At. The third order correction term is 



^(Atf[L A + 2L B ,[L Al L B }}f(Q), 



(40) 



which can be evaluated to estimate the accuracy of the 
approximation. 

In the present case, the evolution operator is L c i + L q + 
Ddp for the Wigner evolution and is the same, but with 

L q = for the dual classical evolution. We split this 
into three operators, the "stream" operator —(p/m)d q , 
the "kick" operator proportional to potential derivatives, 
which differs for the classical and quantum cases, and the 



B. Duffing Oscillator 

The particular potential chosen for study was the 
chaotic Duffing oscillator with unit mass: H(q,p,t) = 
p 2 /2 + Bx 4 — Ax 2 + Ax cos(uit). The evolution was eval- 
uated for the set of parameters A = A = 10, B = 0.5 
and u) = 6.07. In this parameter regime, the system is 
strongly chaotic, with an average Lyapunov exponent of 
A = 0.57 that is relatively uniform over the hyperbolic 
phase space region 4l|. The size of the bounded phase 
space region, which is A in our calculations, is approxi- 
mately 270 units of action. The hyperbolic region of the 
system's bounded motion is generated by the homoclinic 
tangle of a single hyperbolic fixed point and the stable 
regions are relatively small. Consequently the unstable 
manifold associated with the hyperbolic point completely 
characterizes the chaotic region and provides an ideal test 
for the theory developed in this paper for bounded hy- 
perbolic regions. 

These parameters were chosen, not only because they 
provide appropriate testing conditions for theory, but 
also because their classical dynamics have been well stud- 
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ied in Ref. [4J| and elsewhere. As a result, one can esti- 
mate the values of the quantities of interest, e.g., t qc , with 
the system parameters, such as A, fixed at some canonical 
values. In addition, one also must be careful to choose a 
value of h which is not so large that the initial conditions 
are well outside the bounded region. Of course, choos- 
ing a value of h of the same order of magnitude as the 
bounded region or greater, would also invalidate the argu- 
ment. One also does not want to choose D values which 
are very large compared to those at which the transition 
is predicted to occur, as extreme D values, while induc- 
ing quantum-classical correspondence, may wash out any 
intrinsic system dynamics. 

We will principally focus on the case where Ti = 0.1 for 
a variety of practical reasons. The value of H = 0.1 turns 
out to be convenient for these purposes: the critical D 
value predicted is small, but not too small that it is below 
computational resolution, and it also allows a wide range 
of D values to be studied without smearing out the sys- 
tem dynamics. This value of Ti was used in Ref. [l5| which 
motivated much of this research, and which confirms that 
a weak transition will occur for this value. Still, the addi- 
tional set of h values, {0.01, 0.5, 1, V% 3, 5, 10, 20}, were 
studied, and all revealed similar results, though some, 
such as h = 0.01 and Ti = 5, had compromised dynami- 
cal ranges, while % = 10 and 20, were too large to be of 
practical interest. The results presented in depth in this 
section for % = 0.1 should, therefore, be thought of as 
emblematic of all cases studied. For a given trial, h and 
D were held fixed. 

We use the same normalized initial conditions for both 
dual classical and quantum evolutions, as we are trying 
to see the degree to which the two evolutions follow each 
other. In the numerical simulations, the typical condition 
was a superposition of two Gaussians, since a classically 
unacceptable initial condition would better illustrate the 
suppression of interference effects. Other conditions were 
also tried and compared, with analogous results. 

The choice of D = 0.001, yields the estimates t* = 
15.02 and t qc — 57. For this case, the development of 
large scale classical structures should terminate before 
quantum and classical agreement occurs. A larger value 
for the diffusion coefficient, D = 0.01, gives t* = 13.1 
and tqc = 5.7. Here the transition occurs just before 
the termination of classical structure. Because of the 
quantum nature of the initial condition, in the estimation 
of t* , we have set = ft. 

We now set out to test these predictions via numerical 
simulations. We first compare the classical and quantum 
evolutions at late times in order to establish whether or 
not a quantum-classical transition in fact occurs as pre- 
dicted. Comparison of expectation values was helpful to 
establish whether the transition had occurred, but final 
approval was given only after examining the distribution 
functions and Wigner functions directly. Such a compar- 
ison is presented in Figure 2. We compare cross-sectional 
slices of the classical distribution function and quantum 
Wigner function after 149 drive periods of the Duffing 



oscillator. These slices are taken along the p = line. 
For D = 10~ 5 very little agreement occurs between the 
classical and quantum slices (as expected, since in this 
case tqc ~ 5700). In fact, the quantum slice still has 
many negative regions. For D = 10~ 3 , in agreement 
with our order of magnitude estimate for when a transi- 
tion should occur, progress has clearly been made. The 
two functions are in average agreement with one another, 
and, although there is less agreement on details, there is 
agreement between the two on some of the larger phase 
space feature. At a larger value, D = 10~ 2 , this agree- 
ment is much improved and overall, the distributions are 
much smoother. 

Note that the weak QCT occurs in time when l q — 
Id, ~ Vh, this phase space scale being independent of the 
value of D. However, the form of the large-scale clas- 
sical template is determined by the structures present 
at t* , which is sensitive to the value of D. Addition- 
ally, larger values of D will lead to stronger filtering in 
both the quantum and classical cases. At a fixed value of 
time and ti 1 this means that slices of the Wigner function 
at higher D values will have broader features and more 
efficient filtering of small scales. This aspect is clearly 
demonstrated in the three panels of Figure 2. 
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FIG. 2: Sectional cuts of Wigner functions (dashed lines) and 
classical distributions (solid lines) for a driven Duffing oscilla- 
tor, after 149 drive periods, taken at p = for (a) D — 10 -5 ; 
(b) D = 1(T 3 ; (c) D = 1(T 2 . Parameter values are as stated 
in the text; the height is specified in scaled units. 

We now examine the time-dependence in more detail. 
In Figure 3, we display cross-sectional slices taken at 
t = 10 and t = 30. At t = 10, the slice in the top panel, 
the classical and quantum functions have still clearly not 
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explored phase space sufficiently. The Wigner function, 
has significant negative values and the classical distri- 
bution function has not been heavily broken up by the 
dynamics. By t — 30 (still less than t qc ), lower panel, 
the picture begins to resemble the late-time plot shown 
in Figure 2. The negative regions of the Wigner func- 
tion have been largely eliminated and the functions are 
distributed throughout phase space and are in approxi- 
mate agreement. This result is also consistent with the 
estimated value of t nr . 




-0.002 



FIG. 3: Cross-sectional slices of the Wigner function (red) and 
classical distribution function (blue) taken in phase space for 
p = and D — 10 -3 . The higher plot is taken at t = 10 and 
the lower plot is taken at t = 30. 

We now perform a similar analysis for D = 10~ 2 , for 
which t qc is an order of magnitude shorter. The top panel 
of Figure 4 is a snapshot at t = 8, close to t qc , whereas 
the lower plot is taken at the later time t — 20. The 
early-time panel shows that the two distributions closely 
agree on general features of the dynamics, as predicted. 
By t = 20, the weak QCT is well-stabilized, and one sees 
the strong agreement on most individual features present 
in the late-time case shown in Figure 2. As indicated ear- 
lier, similar results were seen at other values of h, at the 
same level of detail shown here for fi = 0.1. Many of 
the studies of the late-time dynamics of this system ap- 
pear in Ref. [23[. As a further example, we show plots 
in Figure 5 for the case of h — 1 and with stronger noise 
coupling than the typical case considered in the theoreti- 




FIG. 4: Cross-sectional slices of the Wigner function (red) and 
classical distribution function (blue) taken in phase space for 
p = and D = 10~ 2 . The higher plot is taken at t = 8 and 
the lower plot is taken at t = 20. 



cal analysis. In the upper panel D = 0.1, while the lower 
plot has D = 1; the snapshots are taken at t = 20. (In 
both cases, t qc <t* .) The top slice shows general agree- 
ment between the two distributions, with some residual 
quantum interference effects. When D = 1, the formal 
value of t qc is less than the dynamical timescale. This 
indicates that the onset of the weak QCT should be very 
rapid and by the relatively late time at which the snap- 
shot is taken, the transition should be complete. The nu- 
merical results are very consistent with this prediction. 
These two plots show that even in the regime where our 
formal analysis might break down (large h, large D), the 
general features and timescales follow the predicted es- 
timates. At even larger (unphysical) values of fi and D, 
the phase space boundaries become important and the 
smoothing gets so large that dynamical features hardly 
survive in the distributions. 

Finally, we address one last point of the argument - 
that the noise-averaged termination of fine scale struc- 
ture would lead to the presence of the early time fold- 
ing associated with the foliation of the unstable manifold 
generated by the homoclinic point of the Duffing sys- 
tem. Evidence for this is presented in Figure 7, a full 
late-time, high-resolution phase space rendering of the 
Wigner function for D = 10~ 3 . The time is taken to 
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FIG. 5: Cross-sectional slices of the Wigner function (red) 
and classical distribution function (blue) following Figure 4. 
The upper panel has D — 0.1 and the lower panel, D — 1. 
Both snapshots are taken for at t = 20. 



be t — 149 (roughly a factor of two greater than t qc ), 
well after the quantum-classical transition has occurred. 
Superimposed on the Wigner function is the early time 
unstable manifold. It is clear that the evolution has or- 
ganized along these early-time features, as expected from 
our analysis. The final distribution which both the dual 
classical distribution and Wigner function approach, once 
the transition has occurred, shows the suppression of the 
late-time, fine-scale features of the unstable manifold, as 
it is supported by the large early-time structures. Quan- 
tum interference, while expected, is local and is strongest 
near the sharp turns in the manifold where branches are 
most close together. This, combined with the previous 
results in this section, allow us to conclude that the basic 
mechanisms posited for the quantum-classical transition 
are consistent with results from numerical simulations. 



VII. CONCLUSION AND FUTURE 
DIRECTIONS 

We have presented a set of physical mechanisms which 
explain the source of the weak quantum-classical transi- 
tion for one-dimensional, bounded chaotic systems. The 
fact that one-dimensional, chaotic systems are being in- 




FIG. 6: Phase space rendering of the Wigner function at time 
t = 149 periods of driving. The early-time part of the unstable 
manifold associated with the noise-free dynamics is shown in 
blue (see text for discussion). The value of D = 10 -3 is not 
sufficient to wipe out all the quantum interference which, as 
expected, is most prominent near sharp turns in the manifold. 



vestigated in real world laboratory experiments, with in- 
teresting potential applications, further enhances the im- 
portance of understanding how unconditioned environ- 
mental interactions affect a subsystem of interest. We 
have used this understanding to derive estimates for the 
time at which the weak transition occurs. 

It is important to keep in mind that currently there is 
no general understanding of which systems will actually 
exhibit a weak QCT, so the existence of this timescale has 
another useful feature, in terms of classification of quan- 
tum dynamical systems. If a weak QCT has not occurred 
by the predicted t qc , our analysis would argue that it 
will not occur at all (within the parametric assumptions 
made). So the existence or nonexistence of this time can 
be used as a test for the occurrence of long-time quantum- 
classical correspondence (but still on timescales shorter 
than the physical equilibration timescale) without any 
knowledge of initial conditions. 

Our numerical results illustrated the compact mani- 
fold structure induced by the bounded phase space re- 
gion. The role of boundedness is a key component in 
the theoretical analysis presented earlier in this paper. 
This topological feature causes the system to fold on it- 
self, which, in turn, allowed us to estimate a timescale 
for the termination of fine structure. A one-dimensional 
bounded chaotic evolution, coupled with noise, appears 
to necessarily terminate fine scale structure. In order for 
the analysis to be valid, the system must be bounded or, 
if it is unbounded, it must at least fold onto itself in such 
a way as to allow a similar process to take place. The 
lack of such an evolution may be a reason why no such 
transition was found for the manifestly unbounded delta- 
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kicked rotor studied in Ref. 18]. Additionally, while an 
unconditioned evolution eliminated all negativity from 
the quantum Wigner function, the distribution was still 
not classical [||. 

It is useful to restate the ways in which the present 
analysis differs from previous work. First, the connec- 
tion to continuous measurement and the weak and strong 
forms of quantum-classical correspondence are explic- 
itly stated. Second, we use the dual classical Fokker- 
Planck equation not to represent a physical classical evo- 
lution, but rather as a dynamical foil of the open-system 
quantum evolution, one that can handle quantum initial 
states, and quantum backaction (which is missing in clas- 
sical theory), but keeps only the classical system prop- 
agator. Third, our analysis is symmetric - we consider 
the effect of noise acting as a filter on the open-system 
quantum evolution (treated in semiclassical approxima- 
tion) melded with a consideration of noise-induced filter- 
ing on the classical dual evolution with its exponential- 
in-time folding of phase space structures characteristic of 
chaos. This folding points to the role of global phase- 
space topology in deriving our results, and distinguishes 
them from local, heuristic analyses of the role of deco- 
herence in the quantum to classical transition . 



Clearly, more work is needed to fully explore the con- 
ditions under which the weak QCT exists, especially the 
role of boundedness. In this regard, investigation of two- 
dimensional systems would be informative as several of 
the topological arguments presented here would likely 
need to be modified. Adding more dimensions would 
introduce effects such as Arnold diffusion which become 
important components of the dynamics. Many of these 
features lack a lower dimensional analog, so it is reason- 
able to believe that they could play an important role in 
the higher dimensional QCT. It would also be interest- 
ing to see what qualitative features of the weak quantum 
to classical transition will be preserved in these systems. 
More immediately, the connection between the require- 
ments for weak and strong QCT scenarios are worth con- 
trasting, especially to delineate parameters regimes for 
validity. This project is currently underway [431 ]. 
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